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Abstract 



In this note, we investigate the regularity of the extremal solution u* for the semilinear 
elliptic equation — Au + c{x) ■ Vw — Xf{u) on a bounded smooth domain of M" with 
Dirichlet boundary condition. Here / is a positive nondecreasing convex function, ex- 
ploding at a finite value a € (0,00). We show that the extremal solution is regular in 
the low dimensional case. In particular, we prove that for the radial case, all extremal 
, solutions are regular in dimension two. 

\^ ' Keywords: singular nonlinearity, advection, extremal solution, regularitjQ. 

in 

, 1. Introduction 

' We consider the elliptic problem 

o ■ 

-Au + c{x) ■ Vw — Xf{u) in fl, 

u>0 inn, (Px) 

^ , u = on dn, 

H , 

Cu ^ where A > 0, f2 is a smooth bounded domain in M" (n > 2), c{x) is a smooth vector field 

over fl and / : [0, a) — > K+ with fixed a e (0, 00) satisfies the following condition (H): 

f is C^, positive, nondecreasing and convex in [0,a) with lim f{t) — 00. 

In the literature, / is refered as a singular nonlinearity. We say that u is a regular solution 
if M G C^(r2), and we also deal with solutions in the following weak sense. 
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Definition 1.1. We say that u is a weak solution of (Px) if < u < a a.e. in n such 
that f{u)d{x,dn) G L^{n) and 

uA(i> - [ Mdiv(0c) = x [ f{u)(f>, V e c^m n H^{n). 

n Jn Jn 

Moreover, u is a weak super- solution of (P\) if " = " is replaced by " > " for all nonneg- 
ative functions (j> G C^(il) n -ffo(^)- 

Clearly, a weak solution is regular if sup^ u < a. For regular solutions, we introduce 
a notion of stability. 

Definition 1.2. A regular solution u of (P\) is said to he stable if the principal eigenvalue 
of the linearized operator L^ x c := —A + c • V — \ f'{u) is nonnegative in H^{fl). 

Exploiting some ideas in [11 ,[13, the solvability of (Fx) is characterized by a parameter 

A*: 

Proposition 1.1. There exists X* G (0, cx)) such that 

• For < A < A*, the problem (Px) has a minimal solution ux, ux is regular and the 
map A I— > liA is increasing. Moreover, ux is the unique stable solution of {Px)- 

• For A = A*, {Px*) admits a unique weak solution u* := lim^-^A* called the 
extremal solution. 

• For A > A*, {Px) admits no weak solution. 

Here the minimal solution means that ux < v for any solution v of (-Pa)- We remark imme- 
diately a close similarity between {Px) and the Emden- Fowler equation with superlinear 
regular nonlinearity, that is 

-Am = \g{u) in C M"; w = on dn, (1.1) 

with A > and g : [0, oo) (0, oo) satisfies 

g is C^, nondecreasing, convex and lim = cx3. (1-2) 

t— foo t 

In fact, there exists also a critical parameter A g (0, oo) for ()l.ip such that all conclusions 
in the above proposition hold true by replacing A* by A (see lul)- It is well known 
by classical examples as g{u) = (1 + u)^ with p > 1 or g{u) = e^, the extremal solution 
u* can be either a regular solution or a real weak solution in the distribution sense with 
sups^ u — oo. 

For general nonlinearity g satisfying (jl.2p . the regularity of the extremal solution u* 
to (II. ip is obtained by Nedev jl3i] for any bounded smooth domain Q, C M" if n = 2, 3; by 
Cabre Q for convex domains in M**; and for radial symmetry case in M" with n < 9 by 
Cabre & Capella {s^. In [l3|, it is proved that, under mild condition on g, the extremal 
solution u* is regular for any smooth bounded domain Q, C M" if n < 9. 
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We can ask the same question about the problem (Pa): For / verifying (H), is it true 
that the extremal solution to (Pa) is regular for general vector field c and general domain 
n C M" with low dimensions n? We will partly answer this question. It is worthy to 
mention that for studying the explosion phenomena in a flow, Berestycki et al. [l| have 
considered the problem (Pa) with a general source / verifying (jl.2[) . 

Without loss of generality, fix a = 1 in the sequel. The problem (Pa) can be linked 
to equation (|l.ip up to the transformation t; = — ln(l — u). In fact, let u solve (Pa), v 
verifies then 

[-Ai>+|Vi'|2+c(a;)-Vw = Ae"/(l-e-^) := Ag(w) in f], 

1 n an ^*5a) 

I V — on OS 2. 

Therefore g verifies (|1.2p and v* = — ln(l — u*) is the extremal solution for the problem 
(Qa)- Thus the regularity of u* is equivalent to the boundedness of v* , however the 
situation could be very different with the presence of advection terms (see 0, [lB|)- In 
last decade, a model describing the steady state of MEMS (Micro-Electro-Mechanical 
Systems) device given by Pelesko and Bernstein in [13] , has drawn many attentions (see 
and the references therein). 

-Au = - — . in C M": u = on dfl. 

(1 — u)^ 

More generally, many precise studies have been done for the singular nonlinearities with 
negative exponent f{u) = (1 — u)^^ {p > 0) in the advection-free situation, i.e. c = 0. In 
that case, when is moreover the unit ball in M", it is known that u* is regular if and 



only if (see [iallO|) 



n<n :^2 + ^+4 /^. (1.3) 

Tending p — 0"*" in (jl.3l) . we see that rip — 2. Therefore we cannot expect in general 
better than dimension two to claim the regularity of u* . 

For the radial case of (Pa), equally when 57 is a ball and c{x) is the gradient of a 
smooth radial function, ux is radial by uniqueness of the minimal solution. We obtain 
the following optimal results which are new even for the advection-free case. 

Theorem 1.1. Assume that n ~ 2, D, ~ Bi. Let j is a smooth radial function and 
c = V7, then the extremal solution u* is regular for any f satisfying {H). 

Theorem 1.2. For any f satisfying {H), D, — Bi and smooth radial function 7, there 
exists C > such that for all A G (0, A*] 

where \ ■ \ is the Euclidean norm in R" . 
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Remark 1.1. The above estimates are optimal. In fact, when f{u) = (1 — u) ^ , p > Q, 
Q = Bi and c = 0, it is well known that u*{x) = 1 — rp+r if n > Up with Up given in 
(11.31). and we have 



n> n„ iff n > 10 or 3 < n < 9, < h 2 + Vn~T. 

But is the extremal solution u* of {P\) regular with general singular nonlinearity 
/ verifying (H), vector field c and smooth bounded domains in R^? The answer is 
affirmative under some additional mild condition on /. 

Theorem 1.3. Assume that f satisfies conditions (H) and the additional conditions, 



(HI) limsup ^ ^ ^ < 1 



and 



fit) 

f'{t){l-t)\n\l-t) 



iH2) limMW^>0. 

Then u* is regular solution to (P\* ) if n — 2, i.e. fl C K^. 

Under more precise conditions on the growth of /, the extremal solution can be showed 
to be regular in some higher dimensions. 

Theorem 1.4. Let f verify (H) and g{v) — e*'/(l — e^"). Assume that g satisfies 
(H3) liminf ^ = l + ,5> 1 

t^oo g{t) 

and 

m liminf^^,>-l-. 

^ ' t^oo g'^{t) 1 + S 

Then v* = — ln(l — u*) is bounded (so u* is regular) when 



1+S 1+S 

Consequently, if fiS > 1, u* is regular for all n < 6. Furthermore, if we can tend 
S to oo, which means g = o{g') near oo, then u* is regular for n < 6 + 'ly/Ji with any 
/X > 0. However, we can never have /i > 1, since otherwise g blows up at finite value 
and contradicts (|1.2p . so the best result we can expect is for n < 9. For example, if 
fiu) — e , then g{v) — e^"*"*^ verifies 5 — co and ji = 1. 

Theorem 1.5. Let f verify (H) and g{v) = e"/(l — e^^). Assume that g = o{g') near 
oo. Rewrite g{t) — g{Q) + ie'''-*-' in (0, oo), suppose there exists to > such that t^h' [t) is 
nondecreasing for t > to, then for any bounded smooth domain D, C M" with n < 9, u* is 
a regular solution. 
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Furthermore, when g = o{g') near oo, the condition {H2) is just equivalent to {H2), 



since 

f"{t)f{t) _{g"-g')g^ ^_ /g"g g \ g 



— ^TT^ ^(s)= ( ) X ( 1-^ ) (s), V^ = l-e-^ 

f'Ht) {g'~gf^' \9'^ 9' J V 9' J ^' 

It is also easy to see that {H3) is equivalent to the condition 

liminf/:Ml^..>0. 

If the equality holds for the whole limit, we have the following optimal result. The case 
f{u) = (1 — u)^^ was obtained in [jj with a different argument. 



u^l- f{u) 



Theorem 1.6. Assume that 



Then u* is a regular solution if n < Up where Up is defined in ([T 

One of the main difhculties here is due to the vector field c{x). When c 7^ 0, the 
operator — A + c- V is not self-adjoint, we use ideas from 7] to get some energy estimates. 
However if c is a gradient, say c = — V7 in il, then — A + c • V can be rewritten as 
where = — div(e'''V) is a self-adjoint operator. In that case, admits a variational 
structure and we can expect more precise estimates of minimal solutions u\, as in the 
radial case. 

The paper is organized as follows: In section 2, we prove quickly Proposition 11.11 and 
show some general consequences of the stability of u\. The section 3 is devoted to the 
proof of Theorems ll.3l to[ 1.6l for general domains. In section 4, we discuss the radial case. 
The norm || • \\q denotes always the standard L'' norm for any q £ [l,oo]. The capital 
letter C denotes a generic positive constant independent of A, it could be changed from 
one line to another. 



2. Preliminaries 

As mentioned above, —A -t- c • V is not a self-adjoint operator for general vector field 
c. However using Lemma 1 in [7|, we have a kind of Hodge decomposition, which tells 
us that for any vector field c G C°°(ri,R"), there exist a smooth scalar function 7 and a 
vector field b e C°°(n,R") such that 

c=-V7 + & and div(e'^6) = in H. (2.1) 

Therefore the problem {P\) can be rewritten as 

- div(e''Vu) + e''b ■ Vu = \e^f{u) in fl. (P() 

On the other hand, we don't have a suitable variational characterization in general to use 
the stability assumption. Fortunately, we can adopt an energy inequality as in Q, which 
is derived from a generalized Hardy inequality of [6|. 
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Proposition 2.2. Let u\ he minimal solution of (Pa)- For any 1 < /3 < 2, we have 

A^eV'(tiA)^'<|^e^|VV.|^ + J^^^eV, ^ e H^{n). (2.2) 
where b is the vector field in (|2.ip . ||6||oo ~ m^-x^ |6(x)|. 

Proof. We use a Hardy type inequality given by Theorem 2 in Q , which says that for a 
positive principal eigenfunction tp of L^^^x.cj for /3 € [1, 2) and any ijj e H^(fl), 



Jn P Jn 

By Cauchy-Schwarz inequality, it is easy to see 



2-/?|V^|2 , 6.Vy^ \b{x)\^ ^ \\b\\\ 



2 (^2 ^ - 2(2-/3) - 2(2-^)' 

so we are done. □ 

Another main ingredient of our approach is just the transformation v = — ln(l — u). 
Let (j> and ^ be nonnegative functions satisfying 0(0) = ^(0) = and = (j)''^. Define 
vx = — ln(l — ux) and g{vx) = e'"^f{l — e~"^). Using {Qx), we get —div{e'^Vvx) + e''b ■ 
Vvx < \e''g{vx) in Q. Let V = (l^ivx) in (E^l), V A S (0, A*), 

A / e'f'{ux^'{vx) 



e 

o 



^Ve(z;A)V«A + Cp / e^02(z;A) 



I / div(eTVf A)e(t'A) + / eT</)2(«A) 



<^ / e'^9ivxmvx)~l I en-avx)^vx + C0 I e'^<l>^{vx) 



/ e''9(vx)avx)+Cp [ e^cj?{vx). 
P Jn Jn 

The last line is due to div(e'''6) ~ 0. We claim then 

Proposition 2.3. Let 1 < /3 < 2. For any A G (0,A*) and any nonnegative test 
functions (j), £, verifying (/)(0) = ^(0) = and ^' — 0'^, there hold 

A / ey'{ux)<P'ivx) < ^ / e'givxmvx) + Cp [ e-'cj.'ivx) (2.3) 



and 



A / ey'{ux)4>\ux)<'^ [ ey{ux)aux) + Cp [ eV(wA). (2.4) 
P Jo Jo 
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The proof of (|2.4p is completely similar to (|2.3p but using (P^) instead of (Qa)- 
We also make use the following behavior of / proved in (l8| . 

Lemma 2.1. i^or any f verifying (H), we have limt^i f(t)/ f'{t) = 0. 

Choose first ^(u) = e" - 1 in (113), then ((u) = and 

Jn P Jo Jo 

Fix 13 e (1,2). By Lemma[2IIl 

A / ey'{ux)e'^^ < C. 
Jo 

Consequently is uniformly bounded, so is ||/(ua)||i- Multiplying (Px) by u\, 

Jo Jo ^ Jo 

which gives 

Proposition 2.4. The family of minimal solutions {ua}o<a<a* is uniformly hounded in 

Remark 2.1. As far as we know, it is always an open question whether the similar 
energy estimation holds for minimal solutions of (11.11) with general regular nonlinearity 
satisfying (jl.2p and general domain VL when n > 6 (see flS,] for n < 5). For the advection- 
free case c — 0, it was proved in [IS] that u* £ H HqITI) under the condition (H), it 
is also true for the gradient case c — V7 (see Lemma \^.l^ . 

Sketches of proof of Proposition [TTTJ We follow the ideas coming from P. [ill. [lo| . 

The main argument is the maximum principle for operators — A + c • V and under the 
Dirichlet boundary condition, we use also the super-sub solution method and monotone 
iteration. 

Let w e Hl{n) be the regular solution of —Ait; + c • Vw = 1 in f2 and fix a > such 
that amaxQ w < \. It is easy to verify that aw is a supersolution of (Pa) for A > small 
enough. As is a subsolution and aw > in fl, (Pa) admits a regular solution for A > 
small enough. As any regular solution u of (Pa) is also a supersolution for (P^) if /i e 
(0, A), the set of A for which (Pa) admits a regular solution is just an interval. Moreover, 
for these A, using [H) and the monotone iteration wq = 0; — Aw„+i + c • Wvn+i — \f(vn) 
in with Vn+i = on d^l for n G N, we get the minimal solution ux = lim„_>.oo w„. 

If we suppose that the principal eigenvalue of L^^^x.c is negative, we can construct, 
as in [J| another solution v < u\ using the associated first eigenfunction, this is just 
impossible by the definition of u\, hence u\ is stable. The uniqueness of stable solution 
comes from Lemmas 2.16 and 2.17 in [8|. 
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Take a positive first eigenfunction of L-^ with tlie Diriclilet boundary condition, by 



A/(0) / e''(p< / Ae^/(u)(p= / \i{L^)uLp ^ / div(e'^6^)u < C. 
Jn Jn Jn Jn 

So A is upper bounded. Define the critical threshold A* as the supermum of A > for 

which (Pa) admits a regular solution, as u* is the monotone limit of u\ when A — > A*, we 

deduce that u* e HQ{fl) is a weak solution of {P\) by Proposition 12.41 

Suppose that -u is a weak solution to {P\). By the monotonicity of /, it is easy to 
verify that for any S > 1, the function v = S~^u is a weak supersolution for {P\/s), then 
the monotone iteration will enable us a weak solution w of {P\/s) satisfying < w < v < 
S^^ < 1. The regularity theory implies then w is a regular solution of {P\/s)- This means 
that X/S < A*. Let 5 tend to 1, we get A < A*. Therefore, no weak solution exists for 
A > A*. 



The uniqueness of the weak solution can be proved in the very similar way as in 11 1 
using the monotonicity and convexity of /, with the strong maximum principle for the 
operator — A + c • V associated to Dirichlet boundary condition, so we omit the details. □ 

3. Regularity of u* for general c and ft 

For proving our results, we will choose suitable functions cj) to apply (|2.3p or (|2.4[) . We 
need also 

Lemma 3.1. For any q > n/2, there exists C > such that the solution v of {Q\) 
satisfies < v < C\\g{v)\\q in fl. 

Indeed, let w be the solution of L{w) :— —Aw + c ■ Vw — Xg{v) in fl with w = on 
dV,. By regularity theory and Sobolev embedding, ||w||oo < C'll^'ll < C' X*\\g{v)\\q 
because q > n/2 > 1. Morover, as L(w — v) > 0, the maximum principle implies then 
< V < w < C\\giv)\\g. 

3.1. Proof of Theorem 1.3 

For simplicity, we omit the index A for u\ or v\. Let 4'{u) — v = — ln(l — u) in p.4|) . so 
^(u) = (1 — u)^^ — 1. Fix /3 G (1, 2) but very close to 2. Repeating the proof of Theorem 



2 in [18| with the assumption {HI), there exists C > such that 

A / e^P^ <C + CCp f e^<l>\u). 



1 

As = o(^(m)) = o(/^) when m 1~, 



Jn 1 - w 



Using the equation {Q\) and d„v < on dil, 

f \Vv\^ = X f e''/(l - e~") + I ^da - f c-Vv<X ( + C\\Vv 

Jn J Jan ov 7^ Jn'^-u 



1 

< C + C||Vw||2. 
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Therefore ||Vu||2 < C, the classical Moser-Trudinger inequality enables us, as n = 2 

/e«"<C„ Vg>l. (3.1) 
Take now (t>{u) — f{u) — /(O) in (|2.4I) . we need to estimate 

Q{u) := f'{u)cty{u) - |e(") = /'(")0(«) - | ^" (s)ds 

= r f'\s)ds-Cf'{u) 



■.= I{u)~-^J{u)~C.f'{u). 

By {H2), there exists 5 > such that 

/(u) - /(O) = r [/'2(s) + r{s)f{s)] ds>{l + S)J{u) - Cf{u), V w e [0, 1) 
Jo 

Let 2^ < /3 < 2, we get C{u) > CI{u) - C. Asserting this in ((^ . 

A / eV'{u)f{u)<C I e'>fiu) + C. 



Consequently, ||/'(u)/^(m)||i < C. By Lemma I^TTl we deduce 11/(^)113 < C. Combining 
with p.ip . ||5(u)||p < C for any p < 3. The proof is completed by Lemma [3.11 as n = 2. 

□ 

3.2. Proof of Theorem 1.4 

Without loss of generality, we can assume that f;(0) — 1. Let (f){t) = g°'{t) — 1 where 
a > is a constant to be determined later. Then 

ft 

12 I 



m = / r (s)rfs 







9^--\s)g'\s)ds (3.2) 



Jo 



2a -I V / 2a - 1 7o 



The condition {H2) yields: Given any e e |^0,/i— xqraj' there exists C > such that 
9{t)g"{t) > (m - - C in [0, oo). Therefore 



*3'"~'(s)5"(5)ds< -(M-e) f g'''-'{s)g'^s)ds + C 
Jo 



(3.3) 



We divide the proof into two cases. 
Case 1: S > 1 and ^ > or 6 < 1 with ^ > ^j^. 
Take a > i. Combine (jS^) and 



2a - 1 



a 



2a -V 



consequently 



2q! — 1 + /i — e 



9 {t)g'{t) + C, for any t > 0. 



(3.4) 



According to (i?3), for any < 6' < S, there exists C > such that g'(t) > {l + S')g{t)~C 
in [0, oo). Setting these estimates in (12. Sp . omitting the index A and recahing that f'{u) = 
g'{v) - g{v), 

5'\ 



T-T7 / eV(«)(5"(«)-ir-CA / e'^{g'^(v)-iy 

< A / e-<nu){g%v) - 1)^ 
2a2A 



< 



;3(2a - 1 + /i - e) 
Consequently, 

5' 



eV"(«)5'(«) + CA / e^ffM + C / 6^(5"^ - 1) 



2^2 



< 



_ 1 + 5' ;3(2a - 1 + /i - e) 
2(5'C 



7 / eV(t')5"(«) + C / e-^g{v) + C / e^(g"(T;) - 1) 



l + (5' 

Choose 5' near (5 such that 

either 5' > 1 and u > or 5' < 5 < 1 with /i > — ^r— • 

1 + 5' ~ A5' 

Through direct computations, for e > sufficiently small and /? = 2 — e, there exists 



a e 



1 6' ^5'{l + 5'){^i-e)-5' 
2' 1 + 5' 1 + 5' 



such that 



For such a, we obtain 



2a' 



1 + 5' I3{2a -1 + ii-e) 



> 0. 



A / eV"M5'W <C, VAe (0,A*). 



(3.5) 



(3.6) 
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Tending now 5' to 5 and e to 0, p.6p holds true provided that 



"<TT^ + — rrs — ■ ^'-'^ 

Therefore 

eV"+^(") < C- / e^9'^''iv)g'{v) + C <C, 

which imphes that ||g(i')||2Q+i < C for a verifying p.7p . Applying Lemma 13. ![ we 
conclude that for n < 2 + 4a with a verifying p.7p . ua is uniformly bounded, hence u* is 
a regular solution if n satisfies (|1.4p . 

Case 2: (5 < 1 and < /i < ii^. 

Now we take a £ — + e), i), the formulas (|3.2I) and (|3.3I) imply then 

The inequality (|3.4p still holds true. Proceeding as for Case 1, we see that for 5' < 5 but 
nearby, e > small and /3 = 2 — e, there exists 



a e , — ^ + ——r, C 



2 '1 + 5' 1 + 5' J V 2 '2^ 

such that (|3.5p is satisfied. Hence we conclude exactly as in Case 1. □ 

3.3. Proof of Theorem 1.5 

Without loss of generality, assume again g{Q) — 1. Take now (t){t) = te"''(*), where 
a > is a constant to be determined, then 

i{t) = f [l + sah'{s)f e'^^'^'Us 

JO 

* [1 + 2sah'{s)] e^^'^'^'Us + f ah^h'^{s)e^"'''^'Us 

Jo 

Thus, for t > to, 
2K{t) 



2a J .s'h'\s)e"'^'^-'Us = C + ^^/^'(s)^ (e^"^-)) 



J to 

where the last integration is considered in the sense of Stieltjes. The monotonicity of s^h' 
in [ioi oo) yields 

K{t) < |t2;j^(^)g2a/i(t) ^ V i > to- 



ll 



So we get 



Using (|2.3p (we drop the index A), 



< 



By Young's inequality, 



1 



/3 

<C I e'' 



s'^l)3/i'(i;)e(2"+l)M") 

Moreover, g = o(g') at infinity yields linit_^oo h'{t) = oo, hence 

^2;j,(^)g2ah(t) ^^3g(2a + l)h(t) ^ 1 1 

i3/t'(i)e(2"+i)'»(t) " 5(<) - 1 ^ 

Fix /3 e {a, 2), the inequality p.Sp implies 

[g{v) - _ f ^,^(2„+i),,(,) < /■ ^,3;,'(„)g(2a + l)h(.) < ^ 



(3.8) 



— >■ as t — >■ (X). 



Recall that g is superlinear, we obtain |l<?(w)||i < C. Consider again w satisfying L{w) 
Xg{v) in and w = on dfl, as w < w in by maximum principle. 



(giv) - 1) 



2q+1 



< c. 



Following the proof of Lemma 2.1 in [17l | (we just need a minor adjustment, say define 
rii = {x E il : g{v) > w^} instead, here T > is a suitable constant), we can obtain that 
if 2a + 1 > n/2, w is uniformly bounded in L°°(n), so does v. Taking 2 > (3 > a > 7/4, 
the result holds for rt < 9. □ 

3.4- Proof of Theorem 1.6 

Here we choose 0(u) = (1 - u)"" - 1 in (HH). For 2A > A* and e > 0, 

2a^ 



P 



- 2e 



-<C, V/3€[l,2). 



I3{2a + 1) "7 ^ (1 - u)P+2«+i 
We have used /'(u)(l - u) > {p ~ t)f{u) - C in [0, 1) by (fTS]) . As e > is arbitrary, 

1 



(1 - u)P+2"+l 



< C 
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provided that 

a2 

■p > , i.e. when a < p + y pip + 1). 

2a + 1 



Therefore \\{l-u)-^\\q < C li q < 1 + 3p + 2^p(p + 1). For any e > 0, as f{u){l-u) < 
{p + e)f{u) + Ce in [0, 1) by ([11]), we have f{u) < C(l - u)-P-\ consequently 

g{v) = eV(l - e-") = < C(l - u)-'-p-\ 

1 — u 

hence |l5(w)|lr < C when 



1 + 3p+2Vp(p+1) 
r < . 

p+l + e 

According to Lemma [STTl the proof is done by taking e — > 0+. □ 
4. Radial case 

As we have mentioned, when c = — V7, the equation (Pa) is rewritten as 

-div(eTVw) = Xeifiu). (4.1) 
With the variational structure, the stability of minimal solutions u\ is equivalent to 

e^lVV^P > A / eV'(«A)^', V V e i/o'(f^)- (4-2) 



Moreover, for any functions and ^ satisfying 0(0) = ^(0) = and ^' = 0'^, the 
estimate (|2.4p is replaced by 

eV'(wA)0'(wA) < / eV(^iA)e(^A). 

Taking now (pit) — f{t) — /(O) and working as for Theorem 1 in [l^, we have 
Lemma 4.1. When c = V7, i/ie extremal solution u* G iJ^ n Hq{Q,). More precisely, 

fMf{ux)<C, VAe(0,A*]. (4.3) 



When 57 = Si is the unit ball, 7(2;) — j{r) with r = ua is radial by uniqueness of the 
minimal solution and satisfies 

r; — 1 

-u" u' --f'u' ^ Xf{u) in (0,1], (4.4) 

r 

with u'{0) = and u(l) = 0. Our main result in this section is the regularity of the 
extremal solution u* for any / satisfying (H) provided n = 2 and the optimal estimate 
for u' claimed in Theorem 11.21 
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The method we use is similar to 0, but the uniform boundedness of ||ma||ci is not 
enough to claim the regularity of it*, because a singular u* could be Lipschitz in many 
cases (see Remark [TTT]) . In fact, the estimate (j4.3|) is crucial for our proof. 

As in since u'^{r) < by maximum principle or equation (j4.4p . the boundedness 

of ||ua|1hi implies that for any fc € N, r > 0, \\ux\\(jk^-g^\g < Ck,r, V A G (0,A*]. So 
we concentrate our attention near the origin. Derivating the equation (14. 4p or (j4.ip with 
respect to r. 



-div (e'^Vu') = e'^u' 



^~ ^ + 'y" 



in (0,1]. 



Using tp — rrj{r)u\{r) as test function in (|4.2p with rj e Hq{Bi) n C{Bi), by similar 
calculation as for Lemma 2.1 in [5j, we obtain 



|V(r7])|^ - (n - l)7f + i'r^r]^ < > 0, V A G (0, A 



(4.5) 



Proof of Theorem 1.1 

For simplicity, we drop the index A. All estimates below hold uniformly for A. First 
as U\ is radial, by maximum principle, we see that u is decreasing in r. Since / and /' 
are nondecreasing functions according to (i?), the estimate (|4.3p implies (as n — 2) 

nr^f'{u{r))f{u{r)) < f f'{u)f{u) <C, V r e (0, 1]. 

By Lemma |2. 11 we have 

C 



f{u{r))<— for ah r e (0,1]. 
r 

Let To e (0, i]. Let be a radial function in Hq{Bi) n C°(i3i) such that 



(4.6) 



ri{r) 



r^^ ifr<ro; 
r^^ ii ro < r < ^, 



and 77 be a fixed function in Bi \i?i/2, independent of tq. The direct calculation yields 



if r < ro; 

if To < r < i. 



Using (|4.5p . as u is uniformly bounded in H^{Bi) by Proposition 
[0, ro] , we get 

' u'{rf 



and r^Tg ^ < 1 in 



Tending tq to 0, there holds 



1 



u'{ry 



-dr < C. 



dr < C. 



(4.7) 
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Consider the following test function used in 15|: For any r < ^ and < < r, 

{(rro)"^ ifs<ro; 
(rs)^^ ii ro < s < r; 
if r < s < i. 

Applying again (|4.5p and combining with ()4.7p . we obtain finally (with vq — > 0) 

^^Li^ds<Cr^, Vr<l. 



(4.8) 



As {e^ru')' = ^\e^rf{u) with n = 2, so e^ru' is nonincreasing in r. Then u'{s) < 
Cru'{r)/s for s e [r, 1], hence m'(s) < Cu'{r) < for any s G [r,2r] if r < i. By 
for any < r < i, 



That means 



2'' „//'„A2 



u'isY 



ds > 



Jr 3 r 



\u{r)\ < Cr in [0,1]. 
However, we need to consider also u" {r) as explained above. Let 

G(r) = e^rw' and *(r) = -2G(Vf) - M / (r - s)/ (M(%/i)) ds 

^0 

where M is a constant to be chosen. Using G' — —Xe^rf{u), 



(4.9) 



Ae^(^)/' {u{s)) ^ + Ae'^(^)/ (zi(s)) - Af/ (zi(s)) 



2s 



< 



< 



2s 

s=-yr 



= 7? 



(V?)) - M / / {u{y/7s)) ds > CA/(0) - CMy^. 



For the last line, we used |7'(s)|/s < C in [0, 1] since 7 is a smooth function (so 7'(0) = 0). 
Fix M > Co + 1, is then concave in [0, 1]. On the other hand, by 



There exists ri > small enough such that vj/' > in [0, ri] with A > Using (|4.4I) . 
(|46)) and (|49)) . for A > ^ and r < ri. 



u (VT") H 1= \-Ju [y/rj 



- CM^/? 



<vI;V)<HM <-2e-(V^)^<C. 
r \/r 
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Applying one more time (14.91) . we see that u"{s/r) > — C for any A > ^ and r < ri. 
Otherwise, by (g^l) and u"{r) < -u'{r)r-^ - i{r)u'{r) < C, wc claim then 

|k"||oo<C, VA>y. 

Combining with (|4.4I) and (14. 9p . it means ||A/(u)||oo < C*, no singularity will occur. □ 

4-. 2. Proof of Theorem 1.2 

As above, we drop the index A and all estimations hold uniformly for A. First, repeat- 
ing the proof of Theorem 1.8, c) in |5], we obtain f'{u{r)) < Cr~^ in (0, 1]. Using Lemma 
[Owith (113]), /(M(r)) < in (0, 1]. Consequently, by (gH), for n > 3, 

nr 

< -e^r"^iM'(r) = / e'^^'h''^^ f{u{s))ds < C s''-^ds < Cr''-^ . 
Jo Jo 

Hence 

|u'(r)| < -. (4.10) 
r 

Let rj he a radial function in Hq{Bi) D C'^{Bi) such that 

r?(r) = 



ro itr<ro; 
^-\/n^ if ro <r <r\. 



in and be a fixed function in B\ \ here rp is any constant in (0, ri), ri > 
is a small constant to be determined. Therefore 



lyy^ m2 / i\ 2 , 2 2 J (7 + 2 - n) Tq ^ ifr<ro; 

(7'V"^ - 2Vn -1 + 1 



2Vn^ 

(7V - 2Vn~T + 1) ifre[ro,ri]. 



We fix ri > small enough such that 



max {7"r^} < min {n — 2, 2^Jn — 1 — l) 

re[0,ri] 

By dl^]), as |V(rry)p - {n - l)rf +j"r^r]^ < for r £ [0,ro], 

r u'\r)r"-^-^^~^dr<C. 

Jro 

Tending rp to 0, we have 

ri 



u'^ry-'-^^^'-^dr < C. (4.11) 



Now we take another test function used in 15|, 

f]{r) = 



if r < ro; 

if tq < r < ri. 
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Combining (|4.5p and (14. lip , we conclude then 



rro 

/ u'2(r)r"-idr < V ro G [0, ri]. 

Jo 

By the monotonicity of e'''r"~^u', similarly as for (|4.9p . it holds 

|u'(r)| < Cr-t+i+v^^ Vre[0,l]. 

Finally, combining with (|4.10p . we are done (in fact, — + 1 + \/n — 1 < —1 for n > 10). 

□ 
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